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In the reply by Ding et alX to our previous comment, 2 
they claimed that the formula we proposed "is a miscon- 
ception resulting from adding an energy and momentum 
dependent function that has to vanish due to fundamen- 
tal reasons" and the formula in their previous work 3-5 
is "valid irrespective of whether the system is infinite or 
extended". In the following, we will demonstrate that 
their claim is wrong due to fundamental reasons. 

They started the discussion from the formula of current 
given by Eq. (2) in Ref.H [see also Eq. (12.21) in Ref.Q, 

I = ^ JdeTr{[f L (s)f L (e)-f R (s)f R (e)][G r (s) 

-G a (e)} + {f L (e)-f R (e)}G<(e)}, (1) 

where all symbols with * represent the matrices in the 
space formed by a complete basic set in the investigated 
system. For graphene, the basic set is usually chosen to 
be the one formed by the joint eigenstates of momentum 
q and pseudospin s (with s — a,b corresponding to A and 
B sublattices, respectively). From the Dyson equation, 
they obtained 

G r (e)-G a (e) = G r (e)[g a {e)- 1 - g r {e)- 1 ]G a {e) 

+ G r (eW(e)-£ a (e)]G a (e). (2) 

Based on the facts that <7°(e) -1 — 5 r (e) _1 = — 2in — > 
and G r,a (e) has neither isolated poles nor a branch cut 
if E r ^ a (e) is finite, they claimed that the first term in 
Eq. ((2J cannot contribute to the current given by Eq. ([T]) 
no matter whether the system is finite or infinite. 

Unfortunately, these authors overlooked an obvious 
fact that all Green's functions and self-energies in the 
infinite system are matrices with infinite dimensions and 
hence the trace over the above term proportional to r\ can 
still be finite. In fact, under the assumption adopted 
in all their works only atoms in the A sublattice are 
connected with the electrodes and the energy and mo- 
mentum dependence of the self-energy can be neglected, 

Le -. S a, Q k s ,k^( £ ) = with a = L,R for the 

left and right leads and N being the number of sites on 
A sublattice, one obtains 

Tr{t a [G r (e)-G a (s)]} 

+ r a g r a (e)(^~j: a )g:(s), (3) 
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FIG. 1: (Color online) |ro<7„(e)| against the energy s for D = 
2.3t 9 (the Debye cutoff) and T\ = Y l L = F = 1.15V 

in which 0J«( E ) = |EQ»( £ ) and 9 r a a (s) = 

qq' 

~k S^qa qa( e )- This means that the contribution to the 
q 

current from the first term in Eq. ^ is proportional to 
the first term in Eq. Q and hence the first term in Eq. (5) 
in Ref. |2j, which has been demonstrated to be nonzero in 
the infinite system in our previous comment^ We would 
like to stress that the above claim is independent of the 
value of the self-energy, in contrast to the claim in the 
reply by Ding et aL— Actually, the first term in Eq. (5) 
in Ref. y is just the spectral density of the system when 
E r,a = 0. However, it is logically incorrect in their reply 
to treat this fact as the evidence that the above term 
becomes zero when E r,a ^ 0. 

We further demonstrate that the first term in Eq @ 
even dominates the current and hence the tunnel magne- 
toresistance (TMR) for parameters used in the work by 
Ding et al.^- From Eq. (5) in Ref. 0, one obtains 

a[ ' a[ ' |l-£T Q (e)l 2 \l-^g r a (e)\ 2 ' 

(4) 

with S r,a = H(Tl +T r )/2. In order to compare the rel- 
ative importance of the two terms in the above equation, 
which just correspond to the two terms in Eq. @, we plot 
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To5a( e )l m Fig- 1 in the case without spin polarization, 
i.e., r£ = = r . Here we adopt the Debye prescrip- 
tion just as Ding et al.&~— the cutoff energy is chosen 
to ensure the conservation of the total number of states 
in the Brillouin zone after linearization of the spectrum 
around the K point, which gives D = 2.3t g . We also 
set To = 0.5D = 1.15t g , which is the largest self-energy 
used in their work*2r— Since this approach and in fact 
any approach with a cutoff are only valid in the energy 
regime far below the cutoff energy D, we only present 
the result with energy \e\ < OAD. From Fig. 1, it is 
seen that |ro3„(e)| is much smaller than 1 in most cases. 
This indicates that the first term in Eq. ((4]) and hence 
the first term in Eq. @ indeed dominates the current and 
hence the TMR. Therefore, due to the erroneous omis- 
sion of the above dominate term, the formula by Ding 
et al£~— is incorrect in the extended infinite graphene 
they investigated and all their results are scientifically 
groundless. In addition, from Eq. one also finds that 

( e ) — Qa ( £ ) ~ 9a ( £ ) ~ 9a ( e ) m this case. This means the 
conductance given by Eq. (18) in Ref. 2 is almost propor- 
tional to the density of states of this system, which is a 
natural result in the weak coupling limit. Therefore, the 
linear dependence on the bias of the conductance from 



our formula, shown in Fig. 1 in Ref. [H, is an adding sup- 
port of the correctness of our result. 

Last but not least, in Fig. 2 in the reply ji Ding et al. 
presented the TMR from the approach using our sug- 
gestion in Ref. d. From their figure, they claimed that 
the TMR from our approach becomes zero at zero bias. 
However, it is impossible to obtain such a result through 
our formula [Eq. (18) in Ref. 0. This can be proved 
as follows. In the zero-bias limit, iTog^e) = 0, thus 
£a(e) " £a(e) = 9 a (s) ~ 9a 0) from Eq. @. Then ex- 
ploiting Eq. (18) in Ref. H and the definition of the TMR 
[1(0) - I(ir)]/I(0) with 1(0) [I(n)} being the current in 
the parallel (antiparallel) configuration of the electrode 
magnetization, one obtains [here7 cr (#) = rjr^(0)/(rj + 
W))] 

E[7<t(0) -7<tW] /pt _ T l\2 
TMR = -2 — = = M: f-. (5) 

E7,(») (rl + rt) 2 

(7 

Obviously, this TMR is nonzero for nonzero spin polar- 
ization (r£ 7^ r^) used in their replyi Apparently these 
authors again made some numerical errors in their com- 
putation. 
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